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ABSTRACT 



We present phase resolved FUV HST FOS spectra of the rapidly rotating, 
highly magnetic white dwarf RE J0317— 853. Using these data, we construct a 
new model for the magnetic field morphology across the stellar surface. From an 
expansion into spherical harmonics, we find the range of magnetic field strengths 
present is 180— 800MG. For the first time we could identify an absorption feature 
present at certain phases at Rill60A as a "forbidden" ls — > 2s component, 
due to the combined presence of an electric and magnetic field. 



Subject headings: White dwarfs — stars: magnetic fields — ultraviolet: stars 
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1. Introduction 

RE J0317— 853 is a highly unusual magnetic DA white dwarf, discovered by Barstow 
et al. (1995, hereafter B95) during the optical follow-up campaign to the ROSAT Wide 
Field Camera (WFC) all sky survey (Pounds et al., 1993). B95 compared the optical 
spectrum and a far-UV (FUV) spectrum obtained with IUE to synthetic spectra generated 
for a number of different field configurations, and found that it was best matched by a 
dipole model, with a polar field strength of «340MG, but which is offset along the dipole 
axis by 0.2 stellar radii in the direction of the southern magnetic pole. In this model, the 
range of field strengths seen by an observer varies between ?«660 MG at the southern pole to 
Ri200 MG at the north pole. A later analysis of circular polarisation data by Ferrario et al. 
(1997, hereafter F97) produced a broadly similar best-fit model: B = 450 MG, displaced by 
35% of the stellar radius and viewed at an angle of 30° — 60° to the dipole axis. 

Fixing the gravity at log g=8.0, the temperature was found, from the shape and 
slope of the continuum in the FUV spectrum, to be 49, 250^^00 K. Thus, not only is RE 
J0317— 853 one of the most magnetic white dwarfs ever found, it is the hottest member of 
this subclass of degenerate objects. 

This unique star set several other records for white dwarfs. The close proximity of an 
ordinary DA white dwarf, LB09802, just 16" away, allowed B95 to estimate the mass of RE 
J0317— 853 at 1.35M , approaching the Chandrasekhar limit. In addition, the discovery 
that it was rapidly and regularly varying in the optical led B95 to conclude that it must 
be rotating with a period of just 725 seconds, making RE J0317— 853 the most rapidly 
rotating isolated white dwarf ever found. Both the fast rotation and the high mass may be 
an indication that RE J0317— 853 is the end product of a double degenerate merger. 

Comparing models to the flux spectrum of a highly magnetic white dwarf can yield 
a variety of good, but not necessarily unique fits (Putney and Jordan, 1995). The optical 
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and FUV flux spectra obtained by B95 for RE J0317— 853 give the field strength over the 
whole star averaged out over several rotations. B95 had no direct information about, or a 
unique solution to, either the direction of the field or the orientation of the magnetic axis. 
Although F97 did attempt to add the circular polarisation spectra in phase, they had no 
knowledge prior to their observations of the optical period, and could not provide a detailed 
model of the field structure across the whole stellar surface. 

In addition, the model described above does not fit the FUV IUE spectrum in detail. 
Here, in contrast to the complicated structure of the optical spectrum where the Balmer 
lines are shifted almost beyond recognition by thousands of A, the mean FUV IUE spectrum 
is comparatively simple. Lyman a is split, with the redward shifted a (lso 2p+i, in 
the field free notation) component being visible in the IUE SWP spectrum. The blueward 
shifted a component lies outside the range of this instrument. However, only one redward 
Lyman a component is predicted by the B95 and F97 models (at 1300 A), yet a second deep 
absorption features is present near 1340A. Such a doublet structure could be explained if 
we are seeing two distinct regions with rather different field strengths. 

This hypothesis could only be tested with the help of phase resolved spectroscopy, by 
which — depending on the angle between the observer and the rotational axis — a major 
part of the global field can be mapped. A 'true' field configuration is one that, when seen 
at different angles due to the rotational phase, matches the data at each observation. In 
this paper we present time-resolved FUV spectra of RE J0317— 853 obtained with the Faint 
Object Spectrograph (FOS) of the Hubble Space Telescope (HST) to look for variations in 
the line positions and strengths of the Lyman a components, and try to determine the field 
morphology across the stellar surface by applying radiative transfer calculations. 
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2. HST Observations and Data Reduction 

FUV observations were obtained across two orbits with FOS on 1996 November 
26, using the blue filter, G130H grating and the 1" aperture, giving approximately lA 
resolution. A total of 142 x35 second exposures were obtained in RAPID mode, 58 on 
the first orbit and 84 on the second. The use of RAPID mode allowed us to minimise the 
readout time between exposures, and maximise the phase coverage across the 725 second 
rotation period. The total exposure time was 4858.6 seconds. The data were extracted 
with standard IRAF/STSDAS routines, and binned into 12 equal phase intervals (Fig. 1). 
Distinct changes in the FUV spectrum are clearly visible across the rotation period. In 
particular, the position of the redward a Lyman a component appears to shift from Ril340A 
at phase 0.167 (phase=0 is defined by our first bin during the HST orbits) to ?«1300A at 
phase 0.667, while the n (lso —> 2po) Lyman a component is also seen to grow in strength 
to a maximum at phase 0.667. One other intriguing absorption feature at «1160A is seen 
to grow in strength at later phases, in conjunction with the £sl300A and stationary Lyman 
a lines. 



3. Analysis and Results 

The observed spectrum and wavelength dependent polarisation of magnetic white 
dwarfs, described by the four Stokes parameters, can be analysed by simulating the 
transport of polarized radiation through a magnetized stellar atmosphere (e.g. Jordan, 
1992, Putney & Jordan, 1995). However, in previous analyses it was extremely time 
consuming to calculate new synthetic spectra for every variation of the magnetic field 
geometry. In particular, when it comes to phase resolved spectroscopy, a different approach 
is necessary, which we have adopted for the first time here. 
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We calculated a library of synthetic spectra for a given effective temperature (40 000 K), 
for different field strengths between 100 and 900 MG (in steps of 1 MG) and for nine 
different angles (equidistant in the cosine) between the magnetic field and the observer. 
Limb darkening was taken into account by a linear approximation based on an interpolation 
between model spectra: = (0.7 + 0.3(j)I((i = 1) (Euchner, 1998). 

In describing the magnetic field geometry we utilised two different approaches: (i) pure 
dipole models with the ability to offset the dipole center relative to the center of the star, 
and (ii) a more general approach involving an expansion into spherical harmonics. 

For a clearly rotating magnetic star, different spectra are observed depending on the 
relative contribution of the magnetic field, and the angles between the observer and the 
magnetic field at a given phase. After prescribing a specific global magnetic field, three 
additional free parameters are needed to calculate a synthetic spectrum for a given phase: 
(i) the angle 90° — % between the rotational axis and the observer, (ii) the angle f3 between 
the rotational axis and the magnetic (z—) axis (relative to which the coordinate system 
of the offset dipoles or the spherical harmonics have been defined) , and (iii) the phase at 
which the magnetic pole passes the meridian between the rotational axis and the center of 
the visible hemisphere. 

The "best" model is found by minimizing the A% 2 between the predicted flux and the 
observations for all the different phases into which the HST observations have been binned. 
The solution is then found by varying the free parameters in a systematic way, using the 
"Amoeba" downhill simplex method in multidimensions (Nelder & Mead, 1965, Press et 
al., 1986). 

In our current analysis we did not vary the effective temperature, although we know 
from the optical (Amy = ±0™1, see B95) and our UV data (Am 1500 ^ = ±0 r . n l) that the 
surface brightness is changing with the rotational phase. This effect may be due to the 
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influence of magnetic pressure on the structure of the atmosphere, but has not been taken 
into account here. Since the fit was slightly better with T cS = 40000 K than with 30000 
or 50000 K we fixed that value and adjusted the theoretical and observed flux to the same 
value at 1500 A. 

Our first attempt was to fit the data with a centered dipole. We were able to roughly 
fit the data with a polar field strength of 220 MG at phase=0.167 or with 450 MG at 
phase=0.667, but not for all phases simultaneously. The next step was the introduction of 
an offset along the magnetic z axis. A reasonable fit was found for a dipole field strength 
of 363 MG, z ff = —0.19 stellar radii, % = 39° and j3 = 20°, very consistent with the solution 
inferred from the (phase averaged) discovery spectrum in B95. 

A slightly better solution was obtained by additional offsets perpendicular to the z 
axis (see Fig. la in Putney & Jordan, 1995): x oS = 0.057, y Q s = —0.04, z oS = —0.220. In 
this case the angles were different, resulting in % — 40° and (3 = 42° (leading to a visible 
range of surface field strengths between 140 and 730 MG); the differences caused by the 
perpendicular offsets were partly compensated by the (5 angle. Since the fit was not perfect 
in either case (reduced x 2 ~ 4 and 3, respectively), we could conclude that there is no 
contradiction to the assumption of a purely axis-symmetrical field. 

This conclusion could also be drawn from an expansion of the surface magnetic fields 

into: 

B r = - E~ 1 E l m =0 (I + 1) (C COS m0 + sin nufy 

P™(cos^) 

Be = + E~ i E l m=0 (9™ cos m0 + hf sin m0) 

dP[ n (cose)/de 
B<p= - E~i EL=o cos + K sin WO 

dP™ (cos 9) /d sin 9 
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with the associated Legendre polynomials P t m . 

For this paper we limited ourselves to I < 3, i.e. to 15 free parameters. The twelve 
parameters with < m < I, and therefore B^, were very small (< 10 MG) compared to the 
axial symmetric contributions from m = (gft = -206 MG , g\ = +77 MG, $ = -39 MG). 
Possibly, slightly better fits can be obtained by including I > 3. However, this does not yet 
seem to be reasonable since we can expect systematic errors introduced by our assumptions 
(e.g. constant effective temperature). 

In this model (with (3 = 29° and % = 34°) the field strengths cover the range 
180-800 MG. However, it is not quite clear whether fields up to 800 MG are indeed present, 
since the the Isq — > 2p +1 component has a maximum (1343 A) at 470 MG, so that higher and 
lower field strengths lead to the same shift. However, both the B95 and F97 analyses of the 
optical and polarimetric data indicate rather strong fields (200-660 MG, and 180-1600 MG, 
respectively). 

Qualitatively, the bimodal behaviour of RE J0317— 853 can be understood in the 
following way: during the "low field" phases (when B < 210 MG dominates, see Fig. 4) 
the ls — > 2p and ls — > 2p +1 components are centered at about 1195 A and 1303 A, 
respectively. When the "high" magnetic fields (B > 280 MG) become more important, the 
ls — * 2p is shifted to smaller wavelengths, out of the HST range, and ls — ► %P+i is close 
to its maximum at 1343 A. 

One unexpected feature, not recognized in the IUE spectrum, is an additional 
absorption line at 1160 A. This appears during the "low field" phase, whenever Isq — > 2po 
is visible as well. In the dense atmosphere of magnetic white dwarf stars, mean electric field 
strengths up to 10 7 V/m can be induced by free electrons and ions. The hydrogen atom can 
experience, due to thermal motion, an additional Lorentzian electric field up to 10 7 • • • 10 9 
V/m, perpendicular to the magnetic field (Mathys, 1989). This gives rise to the motional 
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Stark effect. Without additional electric fields, the magnetic quantum number m and the 
z-parity ix z are conserved. For parallel electric and magnetic fields, only the magnetic 
quantum number m is conserved, and for perpendicular electric and magnetic fields only the 
z-parity. In contrast, no discrete symmetry is left in arbitrarily oriented fields. Therefore, 
in combined electric and magnetic fields, additional allowed dipole transitions occur (e.g. 
Isq — * 2sq, which gains oscillator strength at the expense of the ls — > 2p ) compared to 
the diamagnetic system, where only three transitions (lso - ► 2p-i,o,+i) are dipole allowed. 
Recognizing this possibility, we included the "forbidden" ls — ► 2s in our calculations and, 
consequently, could reproduce the peculiar feature found at 1160 A. 

Our model for the magnetic field is based on the three visible Lyman a components 
(including the "forbidden" one) only. As we have seen, different assumptions about the field 
geometry lead to fits of similar quality. Although the basic properties of the offset dipole 
model and the expansion into spherical harmonics are the same we expect that we will be 
able to constrain the global magnetic field much better when we include, in the near future, 
phase resolved spectropolarimetric data. 

We acknowledge the support of PPARC, UK, DFG (Ko 738/7-1) and DLR (50 OR 
97042). 
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Fig. 1. — Observed and theoretical spectra for 12 different phases across the 725 second 
rotation period. 

Fig. 2. — Distribution of the magnetic field strength over the visible hemisphere from the best 
fitting model for the twelve phases shown in Fig. 1. The steps in grey scale correspond to steps 
of 100 MG, the lowest field is 170 MG (darkest grey), the highest 800 MG (white spot). An 
animation of this plot together with the observed and synthetic spectra can be found on the 
WWW ( |http: / / www.astrophysik.uni- kiel.de/arbeitsgruppen/agkoe/h— jordan_re0317.html|) 

Fig. 3. — Geometry of the best solution model for the expansion into spherical harmonics. 
The upper panel shows the projection of the cone-like movement of the magnetic pole during 
the rotational period onto the celestial plane. The relative orientations of the direction to 
the observer, rotational axis and magnetic pole are shown below. 

Fig. 4. — Relative contribution of low (< 210 MG) and strong (> 280 MG) magnetic fields 
on the visible hemisphere of the model star for the different phases. 
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